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Abstract. We propose a simple and straightforward method based on Wronskians 
for the calculation of bound-state energies and wavcfunctions of one-dimensional 
quantum-mechanical problems. We explicitly discuss the asymptotic behavior of the 
wavefunction and show that the allowed energies make the divergent part vanish. As 
illustrative examples we consider an exactly solvable model and the Gaussian potential 
well. 



1. Introduction 

Wronskians (or Wronskian determinants) are most useful for the analysis of ordinary 
differential equations in generalpQ and the Schrodinger equation in particular [2J. 
Whitton and Connor [3] applied a remarkable Wronskian formalism to resonance 
tunnelling reactions and we have recently proposed that a closely related method may 
be suitable for teaching quantum scattering in advanced undergraduate and graduate 
courses in quantum mechanics[I]. 

The purpose of this paper is to show that the Wronskian method is also useful for 
the study of bound states of one-dimensional quantum-mechanical models. We believe 
that this variant of the approach is also suitable for pedagogical purposes and enables a 
unified treatment of the discrete and continuous spectra of simple quantum-mechanical 
models. 

In section |2j we convert the Schrodinger equation into a dimensionless eigenvalue 
equation and show how to apply the Wronskian method to bound states. In section [3] 
we illustrate the application of the approach by means of two suitable examples. In 
section H] we outline the main results of the paper and draw conclusions. Finally, in 
an Appendix we outline the main properties of Wronskians that are relevant to present 
discussion. 
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2. The Schrodinger equation 

Before solving the Schrodinger equation it is a good practice to convert it into a 
dimensionless eigenvalue equation. In this way one removes all the physical constants 
and reduces the number of model parameters to a minimum. The time-independent 
Schrodinger equation for a particle of mass m that moves in one dimension (—00 < 
X < 00) under the effect of a potential V(X) is 

- |VpO + V(X)^(X) = E^{X) (1) 

where a prime indicates derivative with respect to the coordinate. If we define the 
dimensionless coordinate x = X/L, where L is an appropriate length scale, then we 
obtain the dimensionless eigenvalue equation 

— -ip"(x) + v{x)(f{x) = eip(x) 

<p(x) = VL*J>(Lx), v(x) = ^-V(Lx), e = (2) 

The length unit L that renders both e and v(x) dimensionless is arbitrary and we can 
choose it in such a way that makes the Schrodinger equation simpler. We will see some 
examples in section |3j 

It is well known that a general solution to the second-order differential equation 
(J2J) can be written as a linear combination of two linearly independent solutions. Here 
we write 

<p(x) = A 2 C(x) + B 2 S(x) (3) 

where the solutions C(x) and S(x) satisfy 

C(x Q ) = S'(x ) = 1, C'(z ) = S(x ) = (4) 

at a given point Xq in (—00, 00). These conditions are sufficient to ensure that C(x) and 

S(x) are linearly independent [2] . 

For every value of the dimensionless energy e we know that 

AiL c (x) + BxLdix), x ->■ -00 
<p{x) -> { (5) 

A 3 R c (x) + B 3 R d (x), x ->■ 00 
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where L and R stand for left and right and c and d for convergent and divergent, 
respectively. It means that, for arbitrary e, the wavefunction is a linear combination of 
a convergent and a divergent function when \x\ — >■ oo. If, for a particular value of e, 
B\ = B% = then the resulting wavefunction is square integrable. In other words, this 
condition determines the energies of the discrete spectrum. 

It follows from the Wronskian properties outlined in the Appendix that 

B X W{L C , L d ) = A 2 W(L C , C) + B 2 W(L C , S), x -> -oo 

B 3 W(R C , R d ) = A 3 W(R C , C) + B 2 W(R C , S), x -> oo (6) 

Therefore, when Bi = B 3 = we have a linear homogeneous system of two equations 
with two unknowns: A 2 and B 2 . There will be nontrivial solutions provided that its 
determinant vanishes 

W(L C , C)W(R C , S) - W(R C , C)W(L C , S) = (7) 

The roots of this equation e n , n — 0, 1, . . . are the energies of the bound states (discrete 
spectrum). 

When the potential is parity invariant 

v(— x) = v(x) (8) 

and xo = then C(x) and S(x) are even and odd functions, respectively. In this case 
we have 

W(L CJ S) =W(R C: S) 

W(L C ,C)= -W(R C ,C) (9) 

and the determinant (J7D takes a simpler form: W(R C , C)W(R C , S) = 0. We appreciate 
that the even and odd solutions are clearly separate and their eigenvalues are given by 

W(R C , C) = 

W(R c ,S)=0 (10) 
respectively. Besides, we need to consider only the interval < x < oo. 
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Commonly, it is not difficult to derive approximate expressions for the convergent 
and divergent asymptotic forms of the wavefunction because they are straightforwardly 
determined by the asymptotic behavior of the potential v(x). Therefore, it only remains 
to have sufficiently accurate expressions for C(x) and S(x) and their derivatives in 
order to obtain the eigenvalues by means of equation (|7j). This problem is easily solved 
by means of, for example, a suitable numerical integration method j5]. If y(x) stands 
for either C(x) or S(x) then such an approach gives us its values at a set of points 
%o — N L h, Xq — N L h + h, . . . , x , x + h, . . . , x + N R h where N L and N R are the number 
of steps of size h to the left and right of xq, respectively. The number of steps should be 
sufficiently large so that y(x) reaches its asymptotic value at both xl = xq — N^h and 
x r — x o + N R h and h should be sufficiently small to provide a good representation of 
y(x). The numerical integration methods also yield the derivative of the function y'(x) 
at the same set of points which facilitates the calculation of the Wronskians. 

3. Examples 

In order to test the accuracy of the Wronskian method we first choose the exactly 
solvable problem given by the potential V(X) = — Vo/ cosh 2 (aX), where Vq > and 
a > 0. If we set L = 1/a we are led to the dimensionless Schrodinger equation (121) 
where 



Note that the dimensionless energy e depends on only one independent potential 
parameter vq. The units of length and energy are 1/a and h 2 a 2 /m, respectively, and 
we do not have to bother about the mass of the particle and the Planck constant when 
solving the differential equation. The allowed dimensionless energies are given by[6] 




(11) 




< A - 1 



(12) 
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Figure 1. Wronskians for the exactly solvable problem with e = — 1 and vo = 2.5 

and the spectrum is continuous for all e > 0. It is clear that A — > 1 as vq — > and there 
is only one bound state when 1 < A < 2 (0 < vq < 1). As Do increases more bound 
states appear. As a result there are critical values of the potential parameter for which 
e n = that are given by the condition A n = n + 1 or t> ,n = A n (A n — l)/2 = n{n + l)/2. 

Since limi a .i_ > . 0o i;(a;) = then in this case R c {x) = e~ kx and Rd{x) = e kx , where 
k 2 = — 2e (we only consider the interval < x < oo because the potential is parity 
invariant). Consequently, the allowed energies are determined by the conditions 

W(R C , C) = [C'(x) + kC{x)} e~ kx = 0, x -»• oo 

W(R C , S) = [S'(x) + kS{x)} e~ kx = 0, x -> oo (13) 

for even and odd states, respectively. 

Since the potential (ITTj) is parity invariant we integrate the Schrodinger equation 
from xo = to xr = N^h. Fig. [1] shows the Wronskians W(R C ,C) and W(R C , S) 
for e = — 1 and v = 2.5. We appreciate that x = xr = 5 is large enough to have 
constant asymptotic Wronskians and we choose this coordinate value from now on. 
This numerical test also shows that it is sufficient for present purposes to set h = 0.01 
and Nr = 500 in the fourth-order Runge-Kutta method [5] built in the computer algebra 



system Derive (http://www.chartwellyorke com/derive. html). 



Fig. [2] shows both W(R C , C) and W(R C , S) for x = 5 and v = 6 as functions of 
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Figure 2. Wronskians for the exactly solvable problem with vo = 6 as functions of e. 
Squares mark the exact eigenvalues 
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Figure 3. Wronskians for the exactly solvable problem with e = as functions of vq. 
Squares mark the exact critical parameters 

e. We see that the Wronskians vanish at the exact eigenvalues given by equation f[T2j) . 
This is a confirmation of our earlier assumption that the number of steps and their size 
are suitable for obtaining reasonable results. 

Fig. [3] shows W(R C , C) and W(R C , S) for x = 5 and e = as functions of Vq. In 
this case the Wronskians vanish at the exact critical values Vo jTl = n(n + l)/2. 

Those results for the exactly solvable problem suggest that the Wronskian method is 
successful for all bound-state energies and all well depths. As a second example, in what 
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Figure 4. Wronskians for the Gaussian well with e = as functions of vo 

follows we try the nontrivial problem provided by the Gaussian well V(x) = —Voe~ aX2 , 
where Vo > 0, and a > 0, that we easily convert into the dimensionless potential 

v(x)= - v e~ x \ 

mVo mE 

= T2— , e = -J- (14) 
n a n a 

by means of the length unit L = 1/ y/a. This potential is also parity invariant and 
vanishes asymptotically as \x\ — > oo so that the calculation is similar to the preceding 
example. In order to show that the Wronskian method also applies successfully to 
this model we first obtain some critical values of the potential parameter vq. Fig. H] 
shows W(R C ,C) and W(R C ,S) for x = 5 and e = as functions of vq. The three 
zeros of the Wronskians shown in the figure appear at u 0| i = 1-342, w ,2 = 4.325 and 
^o,3 = 8.898. The numerical methods, although as accurate as one may desire (or the 
computer allows) , are not exact; therefore we can make the coefficient of the divergent 
part of the wavefunction very small but never exactly equal to zero. To illustrate this 
point we consider the ground state of the Gaussian well with vq = 5 and obtain the 
approximate energy e = —3.6077 by means of a naive bracketing algorithm. The 
Wronskian method gives us the coefficient of the divergent part of the wavefunction 
B 3 = W(R c ,C)/W(R c ,R d ) = W(R c ,C)/(2k) = 1.886 x 1(T 6 as a by-product of the 
calculation at x = 5. Note that in this case A 2 = 1 and B 2 = so that C(x) is 
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Figure 5. Unnormalizcd ground-state wavefunction for the Gaussian well with vq = 5 

the unnormalized approximation to (p®. Fig. [5] shows that the wavefunction decreases 
as A 3 e~ kx to almost zero and then increases as B 3 e kx . For intermediate values of the 
coordinate the convergent part dominates because A 3 is much larger that B 3 but for 
larger values of x Rd dominates as shown in Fig. [5j The points in this figure are the 
values of B 3 Rd = 1.886 x io _6 e 2 - 686x . We can improve the calculation and make B 3 as 
small as desired, but it will never be exactly zero because the numerical calculation is 
approximate and the spurious divergent part will always be present. However, in order 
to normalize the wavefunction and calculate expectation values we judiciously truncate 
the coordinate interval at a convenient point and obtain finite accurate values for those 
quantities. The Wronskian method is clearly suitable for deciding on the truncation 
point and estimating the error. 

4. Conclusions 

In our opinion the Wronskian method is sufficiently clear and straightforward for 
teaching an advanced undergraduate or graduate course in quantum mechanics. The 
mathematics requires no special background beyond an introductory calculus course. 
Since many available computer softwares offer numerical integration methods the 
programming effort is relatively light. 




X 
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From a purely theoretical point of view the method is suitable for the discussion 
of the convergent and divergent asymptotic behaviors of the wavefunction and for 
illustrating how the allowed bound-state energies make the divergent part vanish leading 
to square-integrable wavef unctions. The students may try other quantum-mechanical 
models, derive the appropriate asymptotic behaviors analytically and then test their 
results by means of a suitable computer program. They can verify that the Wronskians 
already approach constants as the absolute value of the coordinate increases and that 
the wavefunction already looks like a square integrable function when the coefficient 
of the divergent contribution is almost zero. They can even estimate the remnants of 
the asymptotic divergent part because the Wronskian method provides the necessary 
coefficients. 

In addition to it, the Wronskian method is also suitable for quantum scattering^ 
allowing a unified treatment of both the discrete and continuous spectra of the model. 

Finally, we point out that computer algebra systems are remarkable aids for 
the teaching and learning process because they facilitate the algebraic treatment of 
the problem and even offer the possibility of straightforward numerical calculations 
(although they are considerably slower than specialized numerical programs). 

Appendix A. Wronskians 

In order to make this paper sufficiently self-contained in this appendix we outline some 
well known results about the Wronskians that are useful for the study of ordinary 
differential equations in general [1] and also for the treatment of the Schrodinger equation 
in particular (2j |3]. To this end, we consider the ordinary second-order differential 
equation 



If yi and y 2 are two linearly independent solutions to this equation then we have 



L(y) = y"(x) + Q(x)y(x) = 



(A.l) 



yiL(V2) - y%L(yi) 



d_ 

dx 



W{ yi ,y 2 ) = Q 



(A.2) 
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where 

W(yi, y 2 ) = 2/12/2 - 2/22/1 (A.3) 
is the Wronskian (or Wronskian determinant [1] ) . The Wronskian is a skew-symmetric 

W(f, g) = -W{g, f) W(f, f) = (A.4) 
and linear function of its arguments 

W{f 1 + f 2 ,g) = W(f 1 ,g) + W(f 2 ,g) 

W(cf,g) =cW(f,g) (A.5) 

where c is a constant. 

By linear combination of 2/1 (x) and 2/2(2;) we easily obtain two new solutions C(x) 
and S(x) satisfying 

C(x ) = S'(x ) = 1, C'(x ) = S'(zo) = (A.6) 

at a given point x so that W(C, S) = 1 for all x. If we write the general solution to 
equation (lA.lj) as 

y(x) = AC{x) + BS{x) (A.7) 

then 

A = W(y,S), B = W(C,y) (A.8) 

This equation is quite useful for deriving relationships between the coefficients of the 
asymptotic expansions of the wavefunction in different regions of space as shown in 
sections [2] and [3l Additional mathematical properties of the Wronskians are available 

in Powell and Crasemann's book on quantum mechanics[2]. 
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